Combined analysis of microRNA and mRNA profiles provides insights into the pathogenic resistant mechanisms of the lotus rhizome rot.

Lotus rhizome rot caused by Fusarium oxysporum is a common vascular fungal disease in plants that significantly impacts the yield. However, only a few studies have studied the mechanism of Nelumbo nucifera responding to lotus rhizome rot. Here, we investigated the pathogenic genes and miRNAs in lotus rhizome rot to uncover the pathogenic resistant mechanisms by transcriptome and small RNA sequencing of lotus roots after inoculation with Fusarium oxysporum. GO and KEGG functional enrichment analysis showed that differential miRNAs were mostly enriched in starch and sucrose metabolism, biosynthesis of secondary metabolites, glutathione metabolism, brassinosteroid biosynthesis and flavonoid biosynthesis pathways. Twenty-seven upregulated miRNAs, 19 downregulated miRNAs and their target genes were identified. Correlation analysis found that miRNAs negatively regulate target genes, which were also enriched in starch and sucrose metabolism and glutathione metabolism pathways. Their expression was measured by reverse transcription quantitative PCR (qRT-PCR), and the results were consistent with the transcriptome analysis, thus verifying the reliability of transcriptome data. We selected three miRNAs (miRNA858-y, miRNA171-z and a novel miRNA novel-m0005-5p) to test the relationship between miRNAs and their target genes. The activity of the GUS testing assay indicated that miRNA could decrease the GUS activity by inhibiting the expression of their target genes. Collectively, this study provides a comprehensive analysis of transcriptome and small RNA sequencing of lotus root after inoculation with Fusarium oxysporum, and we identified candidate miRNAs and their target genes for breeding strategies of Nelumbo nucifera.

Dynamic transcriptome profiling provides insights into rhizome enlargement in ginger (Zingiber officinale Rosc.).

The rhizome is an economically important part of ginger (Zingiber officinale Rosc.). However, the mechanism of ginger rhizome enlargement remains unclear. In this study, we performed an integrated analysis of the hormone content and transcriptome of ginger at three rhizome enlargement stages: initial enlargement (S1), middle enlargement (S2), and peak enlargement (S3). With rhizome enlargement, the levels of the hormones zeatin (ZT), gibberellic acid (GA), indole acetic acid (IAA), and jasmonic acid (JA) were significantly increased, and this increase was positively correlated with rhizome diameter. Transcriptomic analysis identified a large number of differentially expressed genes (DEGs); the number of DEGs were 2,206 in the transition from S1 to S2, and 1,151 in the transition from S2 to S3. The expression of several genes related to hormone biosynthesis and signalling and cell division or expansion, and transcription factors was significantly altered, which suggests that these genes play essential roles in rhizome enlargement. The results of correlation analysis suggested that the process of ginger rhizome enlargement may be primarily related to the regulation of endogenous cytokinin, GA3, auxin, and JA biosynthesis pathways and signal transduction; GRAS, HB, MYB, MYB122, bZIP60, ARF1, ARF2, E2FB1, and E2FB2, which may regulate the expression of rhizome formation-related genes; and CYC2, CDKB1, CDKB2, EXPA1, and XTH7, which may mediate cell division and expansion. These results provide gene resources and information that will be useful for the molecular breeding in ginger.

Insight into gene expression associated with DNA methylation and small RNA in the rhizome-root system of Moso bamboo.

Moso bamboo (Phyllostachys edulis), typically a monopodial scattering bamboo, is famous for its rapid growth. The rhizome-root system of Moso bamboo plays a crucial role in its clonal growth and spatial distribution. However, few studies have focused on rhizome-root systems. Here we collected LBs, RTs, and RGFNSs, the most important parts of the rhizome-root system, to study the molecular basis of the rapid growth of Moso bamboo due to epigenetic changes, such as DNA modifications and small RNAs. The angle of the shoot apical meristem of LB gradually decreased with increasing distance from the mother plant, and the methylation levels of LB were much higher than those of RT and RGFNS. 24 nt small RNAs and mCHH exhibited similar distribution patterns in transposable elements, suggesting a potential association between these components. The miRNA abundance of LB gradually increased with increasing distance from the mother plant, and a negative correlation was observed between gene expression levels and mCG and mCHG levels in the gene body. This study paves the way for further exploring the effects of epigenetic factors on the physiology of Moso bamboo.

Comparison of two contrasting Leymus chinensis accessions reveals the roles of the cell wall and auxin in rhizome development.

Leymus chinensis, a perennial native forage grass, is widely distributed in the steppes of Inner Mongolia as the dominant species. The main reproductive strategy of this grass is clonal propagation, which occurs via the proliferation of subterranean horizontal stems known as rhizomes. To elucidate the mechanism underlying rhizome development in this grass, we collected 60 accessions of L. chinensis and evaluated their rhizome development. One accession, which we named SR-74 (Strong Rhizomes), had significantly better rhizome development capacity than the accession WR-16 (Weak Rhizomes) in terms of rhizome number, total and primary rhizome length, and number of rhizome seedlings. Rhizome elongation was positively correlated with the number of internodes in the rhizome, which affected plant biomass. Compared to WR-16, SR-74 had higher rhizome tip hardness, higher abundance of transcripts participating in the biosynthesis of cell wall components, and higher levels of the metabolites L-phenylalanine, trans-cinnamic acid, 3-coumaric acid, ferulic acid, and coniferin. These metabolites in the phenylpropanoid biosynthesis pathway are precursors of lignin. In addition, SR-74 rhizomes contained higher amounts of auxin and auxin metabolites, including L-Trp, IPA, IBA, IAA and IAA-Asp, as well as upregulated expression of the auxin biosynthesis and signaling genes YUCCA6, YUCCA8, YUCCA10, YUCCA11, PIN1, PIN2, UGT1, UGT2, UGT4, UGT10, GH3, IAA7, IAA23, and IAA30. We propose a network between auxin signaling and the cell wall underlying rhizome development in L. chinensis.

Comparative transcriptome analysis on drought stress-induced floral formation of Curcuma kwangsiensis.

The rhizomes and tubers of Curcuma kwangsiensis have extensive medicinal value in China. However, the inflorescences of C. kwangsiensis are rarely known in horticulture, because of its low field flowering rate. In order to improve the flowering rate of C. kwangsiensis, we conducted drought stress treatment on the rhizome of C. kwangsiensis. The flowering rate of rhizome was the highest after 4d of drought stress treatment, and the buds on the rhizome could be obviously swell on the 4th day of rehydration culture. In order to identify the genes regulating the flowering time of Curcuma kwangsiensis, comparative transcriptome analysis was performed on the buds on rhizomes before drought stress treatment, 4 d after drought stress treatment and 4 d after rehydration culture. During this process, a total of 20 DEGs controlling flowering time and 23 DEGs involved in ABA synthesis and signal transduction were identified, which might regulate the flowering of C. kwangsiensis under drought stress. Some floral integration factors, such as SOC1 and FTIP, were up-regulated under drought stress for 4 d, indicating that C. kwangsiensis had flowering trend under drought stress. The results of the present study will provide theoretical support for the application of Curcuma kwangsiensis in gardening.

[Comparison of transcriptome of Atractylodes lancea rhizome and exploration of genes for sesquiterpenoid biosynthesis].

This study compared the transcriptome of Atractylodes lancea rhizome at different development stages and explored genes encoding the key enzymes of the sesquiterpenoid biosynthesis pathway. Specifically, Illumina NovaSeq 6000 was employed for sequencing the cDNA libraries of A. lancea rhizome samples at the growth stage(SZ), flowering stage(KH), and harvesting stage(CS), respectively. Finally, a total of 388 201 748 clean reads were obtained, and 16 925, 8 616, and 13 702 differentially expressed genes(DEGs) were identified between SZ and KH, KH and CS, and SZ and CS, separately. Among them, 53 genes were involved in the sesquiterpenoid biosynthesis pathways: 9 encoding 6 enzymes of the mevalonic acid(MVA) pathway, 15 encoding 7 enzymes of the 2-C-methyl-D-erythritol-4-phosphate(MEP) pathway, and 29 of sesquiterpenoid and triterpenoid biosynthesis pathway. Weighted gene co-expression network analysis(WGCNA) yielded 12 genes related to sesquiterpenoid biosynthesis for the SZ, 1 gene for the KH, and 1 gene for CS, and several candidate genes for sesquiterpenoid biosynthesis were discovered based on the co-expression network. This study laid a solid foundation for further research on the sesquiterpenoid biosynthesis pathway, analysis of the regulation mechanism, and mechanism for the accumulation of sesquiterpenoids in A. lancea.

Differential expression involved in starch synthesis pathway genes reveal various starch characteristics of seed and rhizome in lotus (Nelumbo Nucifera).

Lotus (Nelumbo nucifera) is one of the main aquatic vegetables in China. Its seed and rhizome are main edible parts which are rich in starch. Preliminary experiments of starch contents revealed that seed and rhizome expressed great differences in amylose and amylopectin contents. The rhizomes have higher amylopectin content, while the seeds have higher amylose content. In this study, we have estimated 16 varieties of lotus seeds and found that the amylose content of lotus seeds ranged from 30% to 50%, with an average amylose content of 43%, which showed high-amylose content characteristics. Morphological analysis of lotus seed shown that starch rapid accumulated in 20 DAF (day after fertilization) ∼ 26 DAF. Transcriptome of lotus seeds indicated that starch genes played an important role in seed development. Especially in 22 DAF, the genes which controlled amylose synthesis significantly increased, in contrast, the expression of amylopectin-related genes was stable and might limit the synthesis of amylopectin. We further analyzed the expression patterns of 11 key related genes between lotus seeds and rhizomes, and found that the expression of amylose-related genes to were higher in lotus seed, while the expression of genes related to amylopectin synthesis were higher in rhizome. This study provided a comprehensive research of molecular basis for starch in lotus seed and rhizome. Different expression among key genes during starch accumulation might be the principal cause of the differences in starch properties between seed and rhizome. PRACTICAL APPLICATION: Differential expression involved in starch synthesis pathway genes is the main reason for various starch characteristics of seed and rhizome in lotus. High amylose content in lotus seed is a valuable trait for developing functional food.

Genome and Transcriptome Sequencing Analysis of Fusarium commune Provides Insights into the Pathogenic Mechanisms of the Lotus Rhizome Rot.

Fusarium wilt, a vascular wilt caused by F. commune, has been a serious problem for the lotus. Although some F. commune isolate genomes have been sequenced, little is known about the genomic information of the strain that causes Fusarium wilt of aquatic plants. In this study, the genome of F. commune FCN23 isolated from lotuses in China was sequenced using Illumina and PacBio sequencing platforms. The FCN23 genome consisted of 53 scaffolds with a combined size of 46,211,149 bp. According to the reference genome, F. oxysporum f. sp. lycopersici 4287 isolated from tomato, it was finally assembled into 14 putative chromosomes, including 10 core and 4 lineage-specific chromosomes. The genome contains about 3.45% repeats and encodes 14,698 putative protein-coding genes. Among these, 1,038 and 296 proteins were potentially secreted proteins and candidate effector proteins, respectively. Comparative genomic analysis showed that the CAZyme-coding genes and secondary metabolite biosynthesis genes of FCN23 were similar to those of other Ascomycetes. Additionally, the transcriptome of FCN23 during infection of lotus was analyzed and 7,013 differentially expressed genes were identified. Eight putative effectors that were upregulated in the infection stage were cloned. Among them, F23a002499 exhibited strong hypersensitive response after transiently expressed in Nicotiana benthamiana leaves. Our results provide a valuable genetic basis for understanding the molecular mechanism of the interaction between F. commune and aquatic plants. IMPORTANCE Fusarium commune is an important soilborne pathogen with a wide range of hosts and can cause Fusarium wilt of land plants. However, there are few studies on Fusarium wilt of aquatic plants. Lotus rhizome rot mainly caused by F. commune is a devastating disease that causes extensive yield and quality losses in China. Here, we obtained high-quality genomic information of the FCN23 using Illumina NovaSeq and the third-generation sequencing technology PacBio Sequel II. Compared to the reference genome F. oxysporum f. sp. lycopersici strain 4287, it contains 11 core and 3 lineage-specific chromosomes. Many differentially expressed genes associated with pathogenicity were identified by RNA sequencing. The genome and transcriptome sequences of FCN23 will provide important genomic information and insights into the infection mechanisms of F. commune on aquatic plants.

Integrated mRNA and Small RNA Sequencing Reveals a microRNA Regulatory Network Associated with Starch Biosynthesis in Lotus (Nelumbo nucifera Gaertn.) Rhizomes.

Internode starch biosynthesis is one of the most important traits in lotus rhizome because of its relation to crop productivity. Understanding the microRNA (miRNA) and mRNA expression profiles related to lotus internode starch biosynthesis would help develop molecular improvement strategies, but they are not yet well-investigated. To identify genes and miRNAs involved in internode starch biosynthesis, the cDNA and small RNA libraries of Z6-1, Z6-2, and Z6-3 were sequenced, and their expression were further studied. Through combined analyses of transcriptome data and small RNA sequencing data, a complex co-expression regulatory network was constructed, in which 20 miRNAs could modulate starch biosynthesis in different internodes by tuning the expression of 10 target genes. QRT-PCR analysis, transient co-expression experiment and dual luciferase assay comprehensively confirmed that NnumiR396a down-regulated the expression of NnSS2 and ultimately prevents the synthesis of amylopectin, and NnumiR396b down-regulated the expression of NnPGM2 and ultimately prevents the synthesis of total starch. Our results suggest that miRNAs play a critical role in starch biosynthesis in lotus rhizome, and that miRNA-mediated networks could modulate starch biosynthesis in this tissue. These results have provided important insights into the molecular mechanism of starch biosynthesis in developing lotus rhizome.

A nucleotide signature for the identification of Pinelliae Rhizoma (Banxia) and its products.

BACKGROUND: Ensuring the authenticity of raw materials is a key step prior to producing Chinese patent medicines. Pinellia ternata (Thunb.) Breit. is the botanical origin of Pinelliae Rhizoma (Banxia), a traditional Chinese medicine used to treat cough, insomnia, nausea, inflammation, epilepsy, and so on. Unfortunately, authentic Pinelliae Rhizoma is often adulterated by morphologically indistinguishable plant material due to the insufficient regulatory procedures of processed medicinal plant products. Thus, it is important to develop a molecular assay based on species-specific nucleotide signatures and primers to efficiently distinguish authentic Pinelliae Rhizoma from its adulterants. METHODS AND RESULTS: The ITS2 region of 67 Pinelliae Rhizoma and its common adulterants were sequenced. Eight single nucleotide polymorphisms within a 28-43 bp stretch of ITS2 were used to develop six primer pairs to amplify these species-specific regions. We assayed 56 Pinelliae Rhizoma products sold on the Chinese market, including medicinal slices, powder and Chinese patent medicines, which revealed that about 66% of products were adulterated. The most common adulterants were Pinellia pedatisecta (found in 57% of the assayed products), Arisaema erubescens (9%), Typhonium giganteum (2%) and Typhonium flagelliforme (2%). CONCLUSIONS: A severe adulteration condition was revealed in the traditional medicine market. The species-specific nucleotide assays developed in this study can be applied to reliably identify Pinelliae Rhizoma and its adulterants, aiding in the authentication and quality control of processed products on the herbal market.

Identification of QTLs and a putative candidate gene involved in rhizome enlargement of Asian lotus (Nelumbo nucifera).

KEY MESSAGE: QTL mapping studies identified three reliable QTLs of rhizome enlargement in lotus. NnBEL6 located within the confidence interval of the major QTL cqREI-LG2 is a key candidate gene enhancing rhizome enlargement. Lotus (Nelumbo) is perennial aquatic plant with nutritional, pharmacological, and ornamental significance. Rhizome is an underground lotus stem that acts as a storage organ and as a reproductive tissue for asexual production. The enlargement of lotus rhizome is an important adaptive strategy for surviving the cold winter. The aims of this study were to identify quantitative trait loci (QTLs) for rhizome enlargement traits including rhizome enlargement index (REI) and number of enlarged rhizome (NER), and to uncover their associated candidate genes. A high-density genetic linkage map was constructed, consisting of 2935 markers binned from 236,840 SNPs. A total of 14 significant QTLs were detected for REI and NER, which explained 6.7-22.3% of trait variance. Three QTL regions were repeatedly identified in at least 2 years, and a major QTL, designated cqREI-LG2, with a rhizome-enlargement effect and about 20% of the phenotypic contribution was identified across the 3 climatic years. A candidate NnBEL6 gene located within the confidence interval of cqREI-LG2 was considered to be putatively involved in lotus rhizome enlargement. The expression of NnBEL6 was exclusively induced by rhizome swelling. Sequence comparison of NnBEL6 among lotus cultivars revealed a functional Indel site in its promoter that likely initiates the rhizome enlargement process. Transgenic potato assay was used to confirm the role of NnBEL6 in inducing tuberization. The successful identification QTLs and functional validation of NnBEL6 gene reported in this study will enrich our knowledge on the genetic basis of rhizome enlargement in lotus.

The chrysanthemum DEAD-box RNA helicase CmRH56 regulates rhizome outgrowth in response to drought stress.

Plants have evolved complex mechanisms to reprogram growth in response to drought stress. In herbaceous perennial plant species, the rhizome, which is normally an organ for propagation and food storage, can also support plant growth in stressful environments, and allows the plant to perennate and survive stress damage. However, the mechanisms that regulate rhizome growth in perennial herbs during abiotic stresses are unknown. Here, we identified a chrysanthemum (Chrysanthemum morifolium) DEAD-box RNA helicase gene, CmRH56, that is specifically expressed in the rhizome shoot apex. Knock down of CmRH56 transcript levels decreased the number of rhizomes and enhanced drought stress tolerance. We determined that CmRH56 represses the expression of a putative gibberellin (GA) catabolic gene, GA2 oxidase6 (CmGA2ox6). Exogenous GA treatment and silencing of CmGA2ox6 resulted in more rhizomes. These results demonstrate that CmRH56 suppresses rhizome outgrowth under drought stress conditions by blocking GA biosynthesis.

Transcriptome and metabolite profiling to identify genes associated with rhizome lignification and the function of ZoCSE in ginger (Zingiber officinale).

Ginger (Zingiber officinale Roscoe) is an important spice crop in China, and fresh ginger rhizomes are consumed as vegetable in Sichuan and Chongqing. However, tissue lignification accelerates with rhizome maturation, resulting in the loss of edible quality. To understand the molecular mechanisms of texture modification during rhizome development, we investigated lignin accumulation patterns and identified the key genes associated with lignin biosynthesis using gas chromatography-mass spectrometry (GC-MS), liquid chromatography-tandem mass spectrometry (LC-MS/MS) and RNA-sequencing (RNA-Seq). Results showed that the contents of total lignin and its precursors exhibited notable declines with tissue maturation. However, the lignin composition was remarkably modified and syringyl lignin was deposited in mature rhizomes, leading to ginger lignification. Transcriptome analysis displayed 32 lignin biosynthetic genes were dramatically downregulated with rhizome development, including caffeoylshikimate esterase (CSE ), 4-coumarate-CoA ligase , laccase , cinnamoyl-CoA reductase , cinnamyl-alcohol dehydrogenase , peroxidase and caffeic acid 3-O-methyltransferase , indicating that lignin reduction might be attributed to deficiency in intermediates or the downregulation of key biosynthetic enzymes. Furthermore, overexpressing ZoCSE in Nicotiana benthamiana L. enhanced the total lignin content, suggesting its fundamental role in lignin biosynthesis. RNA-Seq also identified candidate lignin production regulators, including hormone-related genes and NAC/MYB transcription factors (ZoNAC1 , ZoNAC4 , ZoMYB14 and ZoMYB17 ). This result provides a molecular basis for lignin accumulation in ginger.

Hormone and carbohydrate metabolism associated genes play important roles in rhizome bud full-year germination of Cephalostachyum pingbianense.

Cephalostachyum pingbianense is the only woody bamboo species that can produce bamboo shoots in four seasons under natural conditions. So far, the regulatory mechanism of shoot bud differentiation and development is unknown. In the present study, indole-3-acetic acid (IAA), zeatin riboside (ZR), gibberellin A3 (GA3 ) and abscisic acid (ABA) contents determination, RNA sequencing and differentially expressed gene analysis were performed on dormant rhizome bud (DR), growing rhizome bud (GR), and germinative bud (GB) in each season. The results showed that the contents of IAA and ZR increased while ABA content decreased, and GA3 content was stable during bud transition from dormancy to germination in each season. Moreover, rhizome bud germination was cooperatively regulated by multiple pathways such as carbohydrate metabolism, hormone signal transduction, cell wall biogenesis, temperature response, and water transport. The inferred hub genes among these candidates were identified by protein-protein interaction network analyses, most of which were involved in hormone and carbohydrate metabolism, such as HK and BGLU4 in spring, IDH and GH3 in winter, GPI and talA/talB in summer and autumn. It is speculated that dynamic phytohormone changes and differential expression of these genes promote the release of rhizome bud dormancy and contribute to the phenological characteristics of full-year shooting. Moreover, the rhizome buds of C. pingbianense may not suffer from ecodormancy in winter. These findings would help accumulate knowledge on shooting mechanisms in woody bamboos and provide a physiological insight into germplasm conservation and forest management of C. pingbianense.

Differences in gene expression and endophytic bacterial diversity in Atractylodes macrocephala Koidz. rhizomes from different growth years.

Atractylodes macrocephala Koidz. (AMK) is widely used in traditional Chinese medicine owing to its pharmacological activity. Here, we aimed to characterize the differentially expressed genes (DEGs) of one- and three-year growth (OYG and TYG) rhizomes of AMK, combined with endophytic bacterial diversity analysis using high-throughput RNA sequencing. A total of 114 572 unigenes were annotated using six public databases. In all, 3570 DEGs revealed a clear difference, of which 936 and 2634 genes were upregulated and downregulated, respectively. The results of KEGG pathway analysis indicated that DEGs corresponding to terpenoid synthesis gene were downregulated in TYG rhizomes. In addition, 414 424 sequences corresponding to the 16S rRNA gene were divided into 1267 operational taxonomic units (OTUs). Moreover, the diversity of endophytic bacteria changed with species in the OYG (773) and TYG (1201) rhizomes at the OTU level, and Proteobacteria, Actinobacteria, and Bacteroidetes were the dominant phyla. A comparison of species differences among different growth years revealed that some species were significantly different, such as Actinomycetes, Variovorax, and Cloacibacterium. Interestingly, the decrease in the function-related metabolism of terpenoids and polyketides was correlated with the low expression of terpene synthesis genes in TYG rhizomes, as assessed using PICRUSt2. These data provide a scientific basis for elucidating the mechanisms underlying metabolite accumulation and endophytic bacterial diversity in relation to the growth years in AMK.

Transcriptome Analysis of Different Sections of Rhizome in Polygonatum sibiricum Red. and Mining Putative Genes Participate in Polysaccharide Biosynthesis.

As a traditional Chinese herb, the rhizomes of Polygonatum sibiricum Red. are rich in various compounds which have plenty of pharmacological applications and biological activities. Among them, Polygonatum sibiricum polysaccharides (PSP) are the main active ingredients and exhibit a broad range of pharmacological. Based on previous researches, identifying genes involved in PSP biosynthesis will help delineate such pathway at the molecular level. In that case, we performed RNA sequencing analysis for two sections of P. sibiricum Red.'s rhizomes significantly different in PSP content. A total of 435,858 unigenes were obtained by assembling transcripts from both sections and 29,548 (6.77%) ones were annotated in all seven public databases. Analyzing count data of RNA-seq, 13,460 differential expression genes (DEGs) between two sections of rhizomes were acquired. After DEGs were mapped to KEGG databases, twelve represented KEGG pathways related to PSP biosynthesis were summed up. And most DEGs were assigned to the pathway of "Starch and sucrose metabolism". Finally, seventeen candidate genes whose expression levels were related to the polysaccharide content, were considered involving PSP biosynthesis in P. sibiricum Red. The present study lays a foundation of researching the molecular mechanisms of PSP biosynthesis.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing.

KEY MESSAGE: Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnPGGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

Association among starch storage, metabolism, related genes and growth of Moso bamboo (Phyllostachys heterocycla) shoots.

BACKGROUND: Both underground rhizomes/buds and above-ground Moso bamboo (Phyllostachys heterocycla) shoots/culms/branches are connected together into a close inter-connecting system in which nutrients are transported and shared among each organ. However, the starch storage and utilization mechanisms during bamboo shoot growth remain unclear. This study aimed to reveal in which organs starch was stored, how carbohydrates were transformed among each organ, and how the expression of key genes was regulated during bamboo shoot growth and developmental stages which should lay a foundation for developing new theoretical techniques for bamboo cultivation. RESULTS: Based on changes of the NSC content, starch metabolism-related enzyme activity and gene expression from S0 to S3, we observed that starch grains were mainly elliptical in shape and proliferated through budding and constriction. Content of both soluble sugar and starch in bamboo shoot peaked at S0, in which the former decreased gradually, and the latter initially decreased and then increased as shoots grew. Starch synthesis-related enzymes (AGPase, GBSS and SBE) and starch hydrolase (α-amylase and ß-amylase) activities exhibited the same dynamic change patterns as those of the starch content. From S0 to S3, the activity of starch synthesis-related enzyme and starch amylase in bamboo rhizome was significantly higher than that in bamboo shoot, while the NSC content in rhizomes was obviously lower than that in bamboo shoots. It was revealed by the comparative transcriptome analysis that the expression of starch synthesis-related enzyme-encoding genes were increased at S0, but reduced thereafter, with almost the same dynamic change tendency as the starch content and metabolism-related enzymes, especially during S0 and S1. It was revealed by the gene interaction analysis that AGPase and SBE were core genes for the starch and sucrose metabolism pathway. CONCLUSIONS: Bamboo shoots were the main organ in which starch was stored, while bamboo rhizome should be mainly functioned as a carbohydrate transportation channel and the second carbohydrate sink. Starch metabolism-related genes were expressed at the transcriptional level during underground growth, but at the post-transcriptional level during above-ground growth. It may be possible to enhance edible bamboo shoot quality for an alternative starch source through genetic engineering.

Transcriptome profiling unveils the mechanism of phenylpropane biosynthesis in rhizome development of Caucasian clover.

Caucasian clover is the only perennial herb of the genus Leguminous clover with underground rhizomes. However, we know very little about its development process and mechanism. Transcriptome studies were conducted on the roots of Caucasian clover without a rhizome (NR) at the young seedling stage and the fully developed rhizome, including the root neck (R1), main root (R2), horizontal root (R3), and rhizome bud (R4), of the tissues in the mature phase. Compared with the rhizome in the mature phase, NR had 893 upregulated differentially expressed genes (DEGs), most of which were enriched in 'phenylpropanoid biosynthesis', 'phenylalanine metabolism', 'DNA replication' and 'biosynthesis of amino acids'. A higher number of transcription factors (AP2/ERF, C2H2 and FAR1) were found in NR. There were highly expressed genes for R4, such as auxin response factor SAUR, galacturonosyltransferase (GAUT), and sucrose synthase (SUS). Phenylpropanoids are very important for the entire process of rhizome development. We drew a cluster heat map of genes related to the phenylpropanoid biosynthesis pathway, in which the largest number of genes belonged to COMT, and most of them were upregulated in R4.

Exploring the biosynthetic pathway of lignin in Acorus tatarinowii Schott using de novo leaf and rhizome transcriptome analysis.

Acorus tatarinowii Schott is a well-known Chinese traditional herb. Lignin is the major biologically active ingredient and exerts a broad range of pharmacological effects: it is an antitumor, antioxidant and bacteriostatic agent, and protects the cardiovascular system. In the present study, the transcriptomes of the leaf and rhizome tissues of A. tatarinowii Schott were obtained using the BGISEQ-500 platform. A total of 141777 unigenes were successfully assembled, of which 76714 were annotated in public databases. Further analysis of the lignin biosynthesis pathway revealed a total of 107 unigenes encoding 8 key enzymes, which were involved in this pathway. Furthermore, the expression of the key genes involved in lignin synthesis in different tissues was identified by quantitative real-time PCR. Analysis of the differentially expressed genes (DEGs) showed that most of the up-regulated unigenes were enriched in rhizome tissues. In addition, 2426 unigenes were annotated to the transcriptome factor (TF) family. Moreover, 16 TFs regulating the same key enzyme (peroxidase) were involved in the lignin synthesis pathway. The alignment of peroxidase amino acid sequences and the analysis of the structural characteristics revealed that the key peroxidase enzyme had well-conserved sequences, spatial structures, and active sites. The present study is the first to provide comprehensive genetic information on A. tatarinowii Schott at the transcriptional level, and will facilitate our understanding of the lignin biosynthesis pathway.